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A CHARACTERIZATION OF THE NORMAL DISTRIBUTION 
USING STATIONARY MAX-STABLE PROCESSES 


SEBASTIAN ENGELKE AND ZAKHAR KABLUCHKO 

Abstract. Consider the max-stable process r]{t) = maxign 
t £ M'*, where {Ui,i G N} are points of the Poisson process with intensity 
u~^du on ( 0 ,oo), Xi, i G N, are independent copies of a random d-variate 
vector X (that are independent of the Poisson process), and k : R'* —> R is 
a function. We show that the process rj is stationary if and only if X has 
multivariate normal distribution and K{t) — k ( 0 ) is the cumulant generating 
function of X. In this case, 77 is a max-stable process introduced by R. L. 
Smith. 


1. Introduction 


The first class of max-stable processes that was used as a flexible model for 
spatially distributed extreme events was introduced by Smith (1990). Let A be a 
random vector having a zero-mean d-variate normal distribution with covariance 
matrix E. Further, let i G N, be independent copies of A. Independently of 
the Ai’s, let {(7i, i € N} be a Poisson point process on (0, 00 ) with intensity u~^du. 
The stochastic process 


( 1 ) 


Ms{t) := max Ui exp • 




t G 


is now commonly termed the Smith process. This process is max-stable and station¬ 
ary, where the former means that the pointwise maximum - max"^j^ MY,,i of n inde¬ 
pendent copies Ms^i,... ,M-s,n of the process Ms has the same finite-dimensional 
distributions as Ms itself, for all n G N. Recall also that a stochastic process 
G is said to be stationary if it has the same finite-dimensional distri¬ 
butions as the shifted process {M(t + h),t G K'^}, for all h G 


Thanks to its simple form and the small number of parameters in low dimen¬ 
sions, the Smith process has become a widely applied model in spatial extreme value 


statistics (c 

e Haan and Pereira 2006 Engelke et al., 2015 

Oesting 2015 Padoan 

et al. 

2010 

Davison and Gholamrezaee 2012 Westra and Sisson| 2011) and it has 

been extended in several directions 

de Haan and Pereira 

2006 |Kabluchko et al. 

to 

0 

0 

0 

Smith and Stephenson 2009 

Robert, 2013). Smith processes also appeared 


in connection with convex hulls of independent and identically distributed sam- 


pies (Eddy and Gale 

1981 

Hooghiemstra and Hiisler 

1996). It is worth noting 

that the class of Brown-Resnick processes introduced in 

Kabluchko et al. 

(2009) 


includes the Smith process as a special case. In fact, instead of the random parabo¬ 
las {Ni,t) — \ in ([^, it is possible to consider independent copies Zi, * G N, 
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Figure 1 . Realization of the one-dimensional Smith process Q . 


of a zero-mean Gaussian process {Z{t),t G with stationary increments and 
variance function = Ve^rZ{t). The Brown-Resnick process (cf., Brown and 

Resnick [iWf| [Kahluchko et al.| [2^ 


( 2 ) 


MBRit) := max Ui exp 
i^N 




t G 


is stationary, max-stable and its distribution depends only on the so-called vari- 
ogram ^{t) = Var(Z(t) — Z{Q)). Choosing Z to be the random linear Gaussian 
function Z{t) = {N,t) which has cr^{t) = {t, Ef), one recovers the Smith process. 

In the case when the matrix S is non-singular, the Smith process Ms defined 
in PI can equivalently be represented as a moving maxima process (cf., Wang and 
Stoev 2010[ for instance) 


( 3 ) 


Ms(t) = 


det(E) 

{2Tr)d/2 


1/2 


max Vi exp 

iGN 




Here, {{Ti,Vi),i G N} is a Poisson point process with intensity dt x v~^dv on 
X (0, oo). In the moving maxima representation ([^, Ms can be interpreted 
as the maximum over many “storms” where the i-th storm has center point 
and strength Vi , and all storms have a common spatial shape given by the d-variate 
normal density with covariance matrix (as opposed to E in representation 0 ). 
Figure[2shows a one-dimensional stationary Smith process. Note in particular that 
because of stationarity the origin does not play an exceptional role. 

In this note we investigate the following question. If we drop the assumption 
of Gaussianity in 0 and replace the normally distributed random vector N by an 
arbitrary d-dimensional vector X, can we find a function k : —>■ K such that the 
stochastic process 


(4) r]{t) := max Uiexp{(Xi,t) — K(t)}, 

iGN 


is max-stable and stationary? Here, {Ui,i G N} is a Poisson point process on 
(0, oo) with intensity u~^du, and Xi, * G N, are independent copies of X which are 
also independent of the Poisson process. In fact, the max-stability of the process rj 
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follows directly by the construction since @ is the de Haan representation ( |de Haan] 
1984) of the process ry. The crucial part of the question above is the stationarity. 


2. Main result 

Theorem 2.1. The process rj defined in @ is stationary on if and only if X 
has a d-variate normal distribution with some mean p G and covariance matrix 
S, and 


( 5 ) 


K{t) = (/i, ^) + 2 + ^(0)- 


Proof. If X is d-variate normal with mean p and covariance matrix S, and K{t) is 
given by ([^, then it is well known that rj is stationary; see for example Theorem 2 


Kabluchko et al. (2009). 

Let T] be stationary. We have to show that X is normal and ([^ holds. Without 
restriction of generality let k(0) = 0 because otherwise we could replace k(<) by 
K{t) — k( 0) without changing the stationarity of ly. Then, ry(0) has unit Frechet 
distribution, i.e., 

P(d(0) < a;) = P(max[/i < x) = exp{—l/cc}, 

iGN 

for a; > 0. The stationarity of the one-dimensional margins of r/ implies that the 
cumulant generating function 

(6) :=logEe<^’‘>, t e 

is finite, and that K(t) = ip{t), for all t G Indeed, 

P(7y(t) < cc) = P ^max 

= exp I. 

Thus, in fact, K{t) = ip{t) for all t G 

For n S N and ti,... ,f„ G the cumulant generating function of the 

n-variate random vector {{X,ti) — (p{ti ),..., (X,f„) — (p{t„)) is given by 


( 7 ) 


( 8 ) 


• ,Un) := logEexp '^{{X,U) - (p{U))i 
li=i 

/ n \ n 

= (ya ^ ufii - ^ u^ipiti), 


for all Ml,..., Wn S M. Therefore, the general stationarity criterion for max-stable 
processes of the form ([^ given in Proposition 6 in Kabluchko et al.| (2009) implies 
that 7y is stationary, if and only if for all h gW^ and all mi, ..., m„ G [0,1] such that 
Sr=i = 1 it holds that 


( 9 ) 


. ,tn ('^1 : ■ * ■ ; ^n) — +/i,.. (^1 5 • ■ • ; '^n) •> 
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or 


, by 0, 


(10) ip I ^ Uiti j - ^ Uip{ti) = (/? I /i + ^ UiU j - ^ Uip(ti + h). 




2=1 


2=1 


2=1 


For arbitrary ti, t 2 7 ^ and S G [0,1] we have by ( 10 ) that 
P ((1 - 5)ti + 5t2) -(1 - l5)V2(il) - 5^{t2) 

= p{{l — 5)ti + 5t2 + ft.) — (1 — 5)p{ti + ft) — 5p{t2 + ft). 

Applying the above relation with eh instead of ft and rearranging the terms, we 
obtain that for every e > 0, 

p ((1 — S)ti + St 2 + eh) — p ((1 — ft)ti + ^^ 2 ) 


^ _ pjh + eh) - p{ti) ^ ^ p{t 2 +eh) - p{t 2 ) 

£ £ 

Note that the function p is infinitely differentiable by its definition ([^. Sending 
e \ 0 on both sides gives 


( 11 ) 


(Vv? ((1 - 5)ti + 812 ), ft) = (1 - 5){\Jp{ti), ft) + 8 {\Ip{t 2 ), ft), 


where Vp{t) S denotes the gradient of at t G Since equation (111 holds 
for all ft G we obtain 

(12) \/p ((1 - S)ti + St 2 ) = (1 - 8 )Vpih) + 6 Wp{t 2 ), 


for all ti,t 2 G and 6 G [0,1]. Thus, any of the d components of the gradient Wp 
is both convex and concave and hence, an affine function. Consequently, there is a 
d X d matrix S G Matd(ffi) and a vector /i G such that 

Vi^(t) = St + /i. 


Observe that the function ^ —>• M, ^(t) = (/i,t) + ^(t, St), t G has the 

same gradient as p. So, the difference between these two functions has gradient 0 
and we conclude that this difference is constant and, in fact, 0 by the assumption 
(^(0) = k( 0) = 0. We obtain that 

K(t) = p{t) = logEexp{(A,t)} = (/r,t) + i(t,St) t G 

It remains to observe that the matrix S must be positive semidefinite because the 
function p, being a cumulant generating function, is convex. Thus, X has a d- 
variate normal distribution with mean vector p and covariance matrix S, which 
completes the proof. □ 


Remark 2.2. For any random vector X it is possible to make the one-dimensional 
distributions of the process r] defined in (|^ stationary by choosing K{t) to be the 
cumulant generating function of X. However, the above proof shows that except 
in the case when X is normal, it is not possible to make the two-dimensional 
distributions of ij stationary. 


Remark 2.3. Similarly to the setting of the present paper one may ask whether in 
the definition of the Brown-Resnick process © we can replace Z by some more 
general process (for example, a Gaussian process with non-stationary increments). 
This question was studied in Kabluchko (2010[). 
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Remark 2.4. Theorem |2.1| provides a characterization of the multivariate normal 
distribution based on max-stable processes. See (Kotz et ah, 2000 page 151) for a 


review of known characterizations of the multivariate normal law. 
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